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It has been recently shown that randomly charged surfaces can exhibit long range electrostatic 
interactions even when they are net neutral. These forces depend on the specific realization of charge 
disorder and thus exhibit sample to sample fluctuations about their mean value. We analyze the 
fluctuations of these forces in the parallel slab configuration and also in the sphere-plane geometry 
via the proximity force approximation. The fluctuations of the normal forces, that have a finite 
mean value, are computed exactly. Surprisingly, we also show that lateral forces are present, despite 
the fact that they have a zero mean, and that their fluctuations have the same scaling behavior 
as the normal force fluctuations. The measurement of these lateral force fluctuations could help 
to characterize the effects of charge disorder in experimental systems, leading to estimates of their 
magnitudes that are complementary to those given by normal force measurements. 



I. INTRODUCTION 

Stability of soft and biological matter in particular is mostly an outcome of the equilibrium between variable range 
Coulomb and long range van der Waals - Casimir interactions that feature in a plethora of contexts [if. Though 
direct measurements of the latter have been announced by various experimental groups the details and the accuracy 
of experiments are sometimes questioned by the experimentalists themselves Q. Over the last few years there have 
been increasing concerns over how to effectively differentiate between the long range Coulomb interactions and the 
long range van der Waals-Casimir interactions in experiments on interactions between metallic bodies in vacuo A 
number of authors have pointed out that disorder effects may significantly affect the forces between surfaces in such 
experiments; possible sources of disorder include the random surface electrostatic potential [|[ connected with the 
so-called patch effect due to the variation of local crystallographic axes of the exposed surface of a clean polycrystalline 
sample f5j as well as effects of disorder in the local dielectric constant @. The direct detection of disorder effects 
in Casimir force experiments, when the force is measured as a function of the intersurface separation of two plates 
or standardly between a plate and a sphere, is difficult as they must be unravelled from the other forces which are 
always present ||. 

Recently it has been proposed that quenched random charge disorder on surfaces as well as in the bulk can lead 
to long range interactions even when the surfaces are net- neutral @, [1] ■ These long range interactions are induced 
by a subtle image charge effect and they could play a significant role in experiments to measure the Casimir force as 
well as in colloidal science in general [9j, where, for instance, random surface charging can occur during preparation 
of surfactant coated surfaces [Ty] ■ Other examples of objects bearing random charge are random polyelectrolytes and 
polyampholytes 11] . While in the latter case the charge distribution could be quenched (i.e., intrinsic to the chain 
assembly during the polymerization process) or annealed (i.e., when monomers have weak acidic or basic groups that 
can charge regulate depending on the pK of the solution), it is not unequivocal to asses the nature of the charge 
disorder distribution in the case of surfactant coated surfaces. 

In this paper we show that net-neutral surfaces experience two types of disorder generated forces that thus show 
pronounced sample-to-sample fluctuations. The first disorder generated force is normal to the interacting surfaces, 
whose features we have already investigated in detail elsewhere Q, and shows a non-zero average and fluctuations 
proportional to the average. The second one, addressed in detail here, is the lateral disorder generated force, acting 
within the plane parallel to that of the interacting surfaces, whose average is zero but nevertheless exhibits sample-to- 
sample fluctuations which can be quite large. In principle measurements of lateral force fluctuations could be useful in 
characterizing and unravelling the effects of quenched charge disorder and thus help the analysis of its role in normal 
force measurements. 

To give an example of how these forces could be measured one could take a small randomly charged slab at some 
distance I from another slab and place its center randomly at some position opposite the larger slab. The same 
could be done rather more effectively with a plane and a sphere as schematically shown in Fig. [TJ Now due to the 
non-homogeneous quenched charge distribution the smaller slab will experience a random lateral force varying from 
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FIG. 1: A schematic top view of a spherical AFM tip (right: side view) with disordered charge distribution above a planar 
substrate with similar charge distribution. Three different realizations of the experiment, i.e. three different lateral positions 
of the tip above the substrate, are shown corresponding to three different samples of force data. Each sample would show a 
different measurement of the normal as well as lateral force with a sample-to-sample variance calculated in the main text. 



sample to sample. Such forces could conceivably be readily measurable in an SFA type set up [13[ used to measure 
shear forces between solid surfaces sliding past each other across aqueous salt solutions [14[ but with interacting 
surfaces bearing disordered charge distribution. The lateral forces measured in distinct experiments varying in regard 
to the exact relative lateral positions between the interacting surfaces will average out to zero but we predict that the 
fluctuations of this lateral force is non-zero and can give information about the magnitude of charge disorder in the 
system. The fluctuations we compute here thus correspond to sample-to-sample fluctuations and stem from different 
sample (experiment) specific relative positions of the interacting surfaces in different experiments. These fluctuations 
are thus distinct from the temporal fluctuations in the measured force due to thermal fluctuations (an example being 
thermal fluctuations of the instantaneous thermal Casimir force as discussed in [IH). 

Most of our computations are for the slab geometry where they can be carried out exactly; however, we show how 
the lateral force fluctuations can be approximately computed also in the case of the sphere-plane geometry shown on 
Fig. [TJ This configuration is an adaptation of the setup standardly assumed to be within the reach of the proximity 
force approximation (PFA) in the case where the charge disorder on the sphere is assumed to be uncorrelated or 
very weakly correlated. Using an alternative calculational method where there is no dielectric discontinuity (i.e., all 
materials used and the intervening space between them have the same dielectric constant), we can compute the lateral 
force fluctuations exactly for the sphere-plane system. The form of the PFA developed here agrees with this exact 
computation in this limit. For the slab geometry we find that the lateral force fluctuations (lateral force variance) 
behave as A/l 2 where A is the area of the smaller slab and I is the slab separation. In the sphere-plane set up, within 
the PFA we find that the lateral force fluctuations behave as R/l, where R the radius of the sphere, and we take the 
limit where R^S> I, with / the closest distance of the sphere to the plane. 

For completeness we give the expression for lateral force fluctuations in the case where the intervening medium is an 
electrolyte described in the weak-coupling Debye-Hiickel approximation 16j as well. In this case the force fluctuations 
are exponentially screened with a screening length given by the Debye length. 

We then turn to the computation of the normal force fluctuations. The method used here is slightly different 
as in normal force fluctuations there is a contribution from image charges whose average is in general non-zero. We 
reproduce the results of for the average normal force using this method and then go on to analyze its fluctuations. 
For the slab geometry with no electrolyte present, we show that the normal force behaves as A/l 2 , while its variance 
also scales as A/l 2 , making both of them comparable. In the sphere-plane set up, we find that the fluctuations of 
the normal force relative to its average value vary as \JTJR~ and thus become increasingly more important as the 
separation is increased. 

II. LATERAL FORCE FLUCTUATIONS 

Consider two parallel infinite slabs separated by a distance I. The slabs whose surface is at z — has a dielectric 
constant and the slab whose surface is at z = I has a dielectic constant e\. We call these slabs S2 and Si respectively. 
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We denote by e m the dielectric constant of the intervening material. Let each slab have a random surface charge 
density p Q (x) = p a (r,z) with zero mean (i.e., the surfaces are net-neutral) and correlation function in the plane of 
the slabs (r. r' G Si, S2), i.e. 

(p a (r, z)pp(r',z')) =8 a pg as 8(z-l a )5{z' -lp)C a {v-r') a, f3 = 1,2, (1) 

and where we define I2 — and l\ — I. In addition we assume that the charge distribution on slab S\ is restricted to a 
finite area A. In the case where the random charge is made up of point charges of signs ±e of surface density n as then 
we may write g as — e 2 n as , and the correlation function C(r — r') has dimensions of inverse length squared meaning 
that its two dimensional Fourier transform is dimensionless. Typically the values of n s for quite pure samples are 
smaller than the bulk disorder variance which has a typical range of between 1CP 11 to 10~ 6 nm 3 (corresponding to 
impurity charge densities of 10 10 to 10 15 e/cm 3 0]). 
The electrostatic energy of the system is given by 

E=±Jdxcf>(x)p(x) (2) 

where p(x) is the total charge density and 0(x) is the electrostatic potential which is given by 

0(x)= / dyG(x ) y)p(y) (3) 



while G(x, y) is the Green's function obeying 

e V-e(x)VG(x,y) = -S(x-y), (4) 

with e(x) the local dielectric function. Upon changing the charge distribution the corresponding change in the energy 
of the system is thus given by 



SE = J dxdy(Sp(x)G(xj)p(y). (5) 
If pi , the charge distribution on the slab Si , is made up of point charges we have 

Pi( x ) = ^( x - x i)> ( 6 ) 

where qi is the charge at the site x^. Now on moving the smaller slab Si by a distance a laterally, that is to say 
normally to the normal between slabs, we find that the new charge distribution is simply given by 

p'i(x)=^2qiS(x-Xi-&). (7) 

This means that we can write 

<5p(x) = <Jpi(x)=-a-Vp 1 (x). (8) 

As the plate Si is moved laterally the self interaction between the charges in both plates is unchanged, thus the energy 
change is only given by the interaction of the charges and image charges in Si with those in S2. We may thus write 

8E=-a- J dr'drdzdz'V r ,pi(r',z')G{r-r';z,z')p 2 (r,z) (9) 

where r' and r are again the two dimensional coordinates in the planes of <Si and S2 respectively and z' and z' are 
the respective coordinates normal to the planes. We thus note that the integration over the coordinate r' is over a 
finite area A, while that over r is unrestricted. The lateral force on plate Si is thus given by 

8E= -a-F (L) . (10) 

As the charges in plates Si and S2 are uncorrelated we find that 

(SE) = -(a-F( L) ) =0, (11) 
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that is the average lateral force is zero. 

The variance of the energy change is given by 

(5E 2 ) =a i a j { J dr' drdzdz' ds' dsdCdC,' V r <pi(r', z')G(r - r'; z, z')p 2 (r, «)V.jpi(s', (')G(s - s'; C, C)p2(s, C)>, (12) 

where the summation is over the in-plane Cartesian components i, j — 1,2. As the charge distributions on the two 
slabs are independent the only nonzero correlations in the above are given by 

(p 2 (r, z)p 2 (s, 0) = g 2s S(z)S(()C 2 (r - s) (13) 

(V r;/9l (r',z')V s , Pl (s',C')) = g ls 6(z' - l)5(C - l^r.y^ir' - s'). (14) 

This then yields 

(SE 2 ) =a l a J g ls g 2 s J dr'drds'dsG(r - r'; 0, l)G(s - s'; 0, l)C 2 (r - s)V r ;V s ^Ci(r' - s'). (15) 

We now write the above in terms of the two dimensional Fourier transforms, with respect to the in plane coordinates, 
G and C of the functions G and C and carry out the integrations over the unrestricted coordinates r and s to find 

(SE 2 ) = a ^92s J dkdqdrW(Mj G{ kfc 2 {k)C x { q ) e^-KMr'-s'), (16) 

where we have used the fact that G(k) and Ci(k) are functions of |k| — k only. Now using the fact that the surface 
charge patch on slab Si is large (of area A) and assuming that the correlations between charges are sufficiently short 
range we may write 

?2 _ Aaiajgi s g 2s f ^ , ^ n __ n n2/ =, n ^ Aa 2 gi s g 2s 



(SE 2 ) = '7^ 9 s; ^ / dkfc 4 /LG(fc;0,0 Gi(fc)G 2 (fc) = * lflJ "" > / dfcfc J G(fc;0,0 2 Gi(fc)C 2 (fc), (17) 
{2tt) z J 4tt J 

where we have used the isotropy of the k integral. From here we deduce that 

(##) = M f' ! '-' / dfcfc 3 G(fc;0,0 2 Gi(fc)G 2 (fc). (18) 
J 4tt J 

The Fourier transform of the Green's function for the parallel slab configuration is easily computed by standard 
methods and is given by 

G ( k . l) = 2e m exp(-fcZ) 

(,,) fce (e m + ei)(e m + e 2 )(l- AiA 2 exp(-2W)) 



where 



giving the general result in the form 



Aq= £«6™ a = lj ^ (2Q) 



/p(L) F (L) )= AS tj g ls g 2s e 2 n f cxp(-2fcQGi(fc)G 2 (fc) 

Ki 3 ' 7reg(e m + e 1 )2( em + £2 )2 J (i _ Al A 2 exp(-2fc0) 2 ' 1 ' 

When the spatial disorder correlations in both slabs are short range such that G Q (r — r') = S(r — r'), we obtain 

(K iL) ^ L) ) = — A 2n< M %f^ ^ 2A - ln(l - AiA 2 ), (22) 
1 3 4 7 re^ 2 (e m +e 1 ) 2 (e ro + e 2 ) 2 A 1 A 2 

which shows that the lateral force fluctuations decay as A/ 1 2 . We may rewrite this result as 



3 i — 4^,2 2 /2 . , , 



(f> ' f> ') = - ; ij r 2 :: / — , — , (23) 



where the function /(ei/e m , e 2 /e m ) follows directly from Eq. (f2"2")l . It is plotted in Fig. [5] as a function of ti/e m 
and e 2 /e m . As can be easily ascertained, the lateral force fluctuations become weaker as e Q /e m tends to infinity (in 
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FIG. 2: Contour plot of the rescaled lateral force fluctuations, /(ei/e m , ez/^m) (Eqs. (|22[1 and (|23p ), between two parallel 
slabs carrying quenched charge disorder as a function of ei/e m and 62 /e„ 



this case one can see that the function f{x,y) decays, for instance, as f(x,x) ~ lnx/a; 4 and f(x,0) ~ lnx/x 2 when 
x — » oo), which corresponds to the case with perfect metallic slabs. On the other hand, when the dielectric constant 
of the intervening medium is decreased, the force fluctuations become more pronounced and eventually diverge for 
£a/e?n (exhibiting a logarithmic divergence, for instance, as f(x,x) ~ — Inx when x — > 0). 
The above result means that statistically the lateral force behaves as 



ff>~^. (24) 

Another interesting point here is that the lateral force fluctuations are also present when there are no dielectric 
discontinuities in the system. Here if we set e 2 — t\ = e m we obtain the result 

(,) (,) = (25) 

This result can be derived in a rather straightforward but illuminating manner that we derive in the Appendix [A] 

In the case where the intervening medium is composed of an electrolyte with dielectric constant e m and with inverse 
screening length m in the Debye-Hiickel approximation we find that the Green's function obeys 

e V • e(x)VG(x, y) - e e(x)* 2 (x)G(x, y) = -<5(x - y), (26) 

where as before e(x) is only a function of z and k(x) is only non-zero (and equal to a constant k) within the medium 
between the two slabs. From this we obtain 

Q (k . q n 2e m Kexp(-Kl) 

{ ' ,} e (e m K + e 1 k){e m K + e 2 k){l-A lK A 2K exp(-2Kl)) [ ' 



where K = \Jk 2 + k 2 and 



A ^ = e a k e m K a = ^ (28) 



In order to obtain the force fluctuations for a system with an intervening electrolyte, at the level of the Debye-Hiickel 
approximation, we simply need to use the expression (|27[) in Eq. ([T| 



A. PFA for lateral forces 



In many experimental set ups, due to problems of achieving a perfectly parallel alignment, a sphere-plane configu- 
ration is used rather than a plane-parallel configuration. In the case where there is no dielectric discontinuity we can 
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compute the lateral force fluctuations for the sphere-plane geometry. The derivation given above is easil y m odified to 
the case of general geometries if one assumes the validity of the proximity force approximation (PFA) 15]. We find 
in that case that for the sphere-plane geometry the force correlator is given by 

4 J 327r 2 e2 e 2j [(x-y)2 + z(x,y)2] 3 ' 

where x are the Cartesian coordinates on surface Si of object 1 (here the surface of a sphere of radius R) projected 
onto the in- plane coordinates of surface 2 and y the coordinates on surface 2 or S 2 (here an infinite plate). The 
variable z(x, y) is the distance between the points on the two surfaces perpendicular to the surface S 2 - In terms of 
spherical polar coordinates on the surface Si if x = (R sin 9 cos ip, R sin 8 sin ip) then we have z(x, y) = I + R(l — cos 8) , 
where I is the distance between the opposing pole of Si and the plane S2 (or the closest distance of the sphere to the 
plane). The integral can now be written as 

(P f = ^-f- f R 2 sintf^dz -3, (30) 

■'-"'''«'■» 1 La + (l + R(l-cos8)Y 



where z is the relative coordinates of Si and 62 in the plane of S 2 (i.e., it represents x — y where x is in the plane of 
Si and y in the plane of S 2 ). Performing the integral over z we then find 

(ff >if >) = 6l f o ff LinedO l - - 2 (31) 

1 J 32e 2 ,e 2 ra J (I + R(1- cos 8)f 

and finally the integral over 8 is easily carried out to give 

<*< *' } ~ 1^1(1 + 2RY (S2) 
In the usual experimental set up we are in the limit where R 3> I and we thus find 

<*f >if >> - (33) 

In the case where there are dielectric discontinuities we can try to approximate the computation of the force correlator 
in a manner similar to the proximity force approximation for electrostatic and Casimir interaction problems. When 
the charge distribution are delta-correlated we can assume that the force due to the interaction of a unit of area on 
the sphere at the same separation from the plane (thus a ring on the sphere) is statistically independent of the others. 
The ring is specified by the polar angle 8 and using Eq. (|22p we can write that the force on a ring of polar angle 
between 8 and 8 + 88 is given by 



pWrm- / 9xs9is^ n a A ) y ^(0)V2ttR 2 sin 8 88 

Fl (0) -r^el(e m + eine m + e 2 rAiA 2 Hl - AlA2)X l + R(l-cos8) ' (34) 

where all the prefactors /i(#) are independent and are of zero mean and variance one. The correlation function of the 
total force is thus given by 

^Ff) = A ^ tr ^g^T^ A _ -(I - Ar A 2 ) / „ (35) 

which gives 



47re2(e m + ei) 2 (e m + e 2 ) 2 AiA 2 'J [I + R(l - cos 8)} 



{Fi F > } ~ ~el(e m + ei)He m + e 2 ) 2 A 1 A 2 ln(1 " ^WTwj ' (36) 



and clearly corresponds to the exact result Eq. (I3"2l in the case where there are no dielectric discontinuities. 
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III. NORMAL FORCE FLUCTUATIONS 

The magnitude of the normal force has been obtained previously 0, [1] and we concentrate our efforts to its 
fluctuations. The calculations for the normal forces between dielectric slabs with random surface charging are slightly 
different to those above for lateral charges. Here we proceed by writing the electrostatic energy as 

E=\J dxdy p(x)G(x, y; l)p(y) (37) 

where we have made explicit the dependence of the Green's function on the slab separation I. The electrostatic 
component of the force of the slabs in the normal direction is then given by 

F (JV) = i [dxdyp(x)H(x,y;l)p(y) where ff(x,y,0 = -^G(x,y;Z). (38) 

The average value of the normal force is non-zero due to the correlation between the charges in each plate and their 
image charges [7|. In terms of the notations introduced earlier we find 



(F {N) ) = -^ J kdk [ 5ls #(fc;0,0)<5i(fc) + g 2s H(k;l,l)C 2 (k)]. 



(39) 



The two terms above are the interaction of the charges on surface 1 and 2 with their images. Note that the contribution 
of the two surfaces are additive as they are independent. 

In the case where there is electrolyte in the region between the two plates that can be described on the Debye - 
Hiickel level one again obtains the relevant expressions for the Green's functions above as 

r(k n a) = 1 ( 1 - A lK cxp(-2^Q \ 

UlK ' U,UJ e (Ke m + ke 1 )\l-A lK A 2K exp(-2Kl)J [W) 
G(k- 1 I) = 1 ( l-A 2K exp(-2jfQ \ 

U[K ' l ' l> e (Ke m + ke 2 )\l-A lK A 2K exp(-2Kl))> [U) 

and this then gives the corresponding derivatives H as 

H(k; 0, 0) = ^e m A 2K e M -2Kl) 

e (Ke m + ket) 2 (1 - A 1k A 2k exp(-2Kl)f 

H(k; l, l) = ^e m A lK e M -2Kl) 

e (Ke m + fce 2 ) 2 (1 - A 1k A 2k exp(-2Kl) f 

H(k; 0, = 2^ em exp(-^)(l + A lK A 2K exp(-2XQ) 

e {Ke m + ke 1 )(Ke m + ke 2 ) (1 - A 1k A 2k exp(-2^)) 

The definition of A aK was given in Eq. 1281 The normal force fluctuations may be computed using Wick's theorem 
and are given by 



(F {N)2 )c kdk [g 2 ls H 2 (k; 0, 0)Cf(k) + g 2 2s H 2 (k; I, l)C 2 (k) + 2g ls g 2s H 2 (k; 0, l)C x {k)C 2 {k)\ . 



(45) 



The first two terms are the force fluctuations due to the self interactions, i.e. of the charges on surfaces 1 and 2 with 
their images, and the last term is the fluctuations of the force between the charges on surface 1 with those on surface 
2 (whose average is always zero). 

If we take the limiting case where the intervening medium is a simple dielectric devoid of any electrolyte, i.e. k = 0, 
and where the surface charges are not spatially correlated so that C a (r — r') = 8(r — r'), we find that the average 
value of the normal force is given by 

w Ae m ln(l-A 1 A 2 ) / g ls g 2s \ 

{ ' 4W 2 vAx^ + ex) 2 + A 2 (e m + e 2 y) ' 1 ' 

which recovers our previous results for the average of the normal force due to quenched charge disorder 0, Q . This 
result may be rewritten as 

{FW) ^Agu gf^^iuy (47) 



47re e m ? 2 \£m Cm gis 
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where the function G(ei/e m ,e2/e m i92s/9is) follows directly from Eq. (|46l) and is shown in Fig. [3] for the case with 
9is — 92s- Note that in this case the average normal force changes sign and turns from repulsive to attractive when 
ci/em and £2/e m become larger than a certain value (shown in the figure by the contour line labeled by 0). For the 
symmetric case with Ai = A2 = A, one has an attractive force when A > (e.g., for two dielectric slabs interacting 
across vacuum) and a repulsive force when A < 0. The normal force diverges logarithmically when e a /e m — > as well 
as when both dielectric constants t\ and e 2 tend to infinity (perfect metal limit). However, it can take a finite value 
when only one of the dielectric constants tends to infinity (note, for instance, that Q(x,0) — > — In 2 ~ —0.69 when 
x — > 00). 

In this case the normal force fluctuations variance are given by 



A 



47Tf 2 e 2 P 



(gf s D n + gl s D 22 



2gi s g2sD 2 i) 



(48) 



where 



D n = 



D22 = 



D21 = 



2ef 



3(e m + ei) 4 
2ei 



3(e„ 



£2 



Ai(l- AiA 2 



Ai 



A 2 (l- A!A 2 ) 2 



3(e m + ei) 2 (e m + 



AiA; 



ln(l- AiA 2 ) 
A? 

ln(l - AiA 2 ) 

A 2 
^2 

■ ln(l - A]A 2 ) 



(I-A1A2 



(49) 
(50) 
(51) 



These expressions for Dn and D21 are shown in Figs. 2£i and b as a function of e\/e m and £2/e m (note that -D22 
can be obtained from D\\ by replacing the subindex 1 with 2 and vice versa). In Fig. |2p, we show the quantity 
Dn + D22 + 2£>2i which can be defined as the rescaled normal force fluctuations for the case with g ls = g 2s through 



A 9u 



Att^ 2 12 



( ni 



£m ' 9li 



(52) 
0, i.e., 



The different contributions D a p to the normal force fluctuations all diverge algebraically when e a /t m 
C(x, x) — > x~ 2 when x —> 0. 

In the case where there are no dielectric discontinuities the forces due to image charges are zero and the only normal 
force is due to the interaction between the charges on the two (net-neutral) surfaces (one can easily see that D u = 
when 62/e m = 1; same is true for D22 — when ei/e m = 1, which explains the non-monotonic behavior of Du and 
£(eiAm, £2/em, 1) as seen in Fig. and c). The mean of the normal force is clearly zero in this case but it has a 



9 




FIG. 4: Contour plots of (a) log 10 Dn, the contribution to the force fluctuations due to the self interactions, i.e. of the charges 
on surface 1 with their images, and (b) log 10 D 2 i , the contribution to the force fluctuations due to the interaction between the 
charges on surface 1 with those on surface 2 as a function of ei/e m and £2 Am- (c) shows the rescaled normal force fluctuations 
log 10 £ (ei/e m , e 2 /e m , g 2a /g al ) for g la = g 2a . All plots are shown in log 10 - log 10 scale for ei/e m and e 2 /e m . 



non-zero variance 

< FW ^3^^' (53) 

a result which can be verified using the expression for the Coulomb potential in a system of constant dielectric constant 
e m , with a computation similar to that leading to Eq. (|A4I) and then Eq. (|25[) . Interestingly we see, comparing with 
Eq. ([25]) . that in the case of a uniform dielectric constant the variance of the force fluctuations in the normal direction 
are twice the magnitude as those in the lateral direction. 



A. PFA for normal forces 



Within the proximity force approximation for the sphere-plane geometry in complete analogy to the case of lateral 
force we can derive both the normal force as well as its fluctuations. The former can be obtained in the form 

{N) = i? 2 e m ln(l-A 1 A 2 ) / g ls g 2s \ 

[ ' e l(l + 2R) {A^ + erf + A 2 (e m + e 2 )^J ' 1 ' 
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and the normal force fluctuations as 



I? 2 



(55) 



elell(l + 2R) 



From this formula we see that the relative size of the fluctuations of the normal force to its average scales as 



and thus the fluctuations of the normal force relative to its average value become more important as the separation is 
increased ! The sample-to-sample scatter in the normal force thus increases on increasing the separation between the 
interacting bodies. This could be interpreted to mean that the interactions themselves restrict their own fluctuations. 



In this work we have proven that the sample-to-sample variance in the lateral as well as normal charge disorder 
generated forces can be substantial. In the ideal, thermodynamic type, limit where the probe area is very large the 
force is much large than its fluctuations. However in some experimental set ups the probe size may be quite small 
and so sample to sample force fluctuations could become important with respect to average forces. In addition we 
have shown that for the sphere plane set up fluctuations become important at large separations where the normal 
force is weak. In the case of lateral force variance, since the average is, the fluctuations are the only thing remaining. 
Interestingly enough the fluctuations in the normal and lateral direction are always comparable. For the special case 
of a uniform dielectric constant we also showed that the variance of the force fluctuations in the normal direction is 
exactly twice the magnitude of the one in the lateral direction. 

The sample-to-sample variation in the disorder generated force is fundamentally different from the thermal force 
fluctuations in (pseudo)Casimir interactions as analyzed by Bartolo et al. [Hj]. In this case one could (in principle at 
least) use the same experimental setup and just observe the temporal variation of the force at a certain position of 
the interacting surfaces, measuring the average and the variance within the same experiment (assuming a sufficiently 
good temporal resolution of the force measuring apparatus). On the other hand, in order to detect sample-to-sample 
variation one would have to perform many experiments and then look at the variation in the measured force between 
them. In the first case the variance of the force is intrinsic to the field fluctuations, in the second one it is intrinsic to 
the material properties of the interacting bodies. 

Additionally, the variance of the fluctuation-induced Casimir force is not universal and is intrinsically related to 
the microscopic physics that governs the interaction between the fluctuating (elastic in the case investigated in [12[) 
field and the bounding surfaces. 

There are several assumptions in our calculation that need to be spelled out explicitly. We always assume that 
the position and orientation of the interacting surfaces are fixed in these experiments as well as in the corresponding 
calculation, just as indicated in the schematic representation of our system on Fig. [TJ However, for an unconstrained 
colloid particle rotational degrees of freedom are not quenched but rather annealed. For example a spherical colloid 
will rotate so as to minimize its interaction energy in the same way as permanent dipoles orientate with each other. 
This would introduce additional considerations in the analysis of forces that we do not address in this contribution. 
In principle there will be random torques and their sample-to-sample variation can be computed using the methods 
presented here. These random torques may be accessible to torsion balance based setups jl7j |. 

Additionally, in force measurements slight translations of the sphere with respect the plane and slight rotations of 
the sphere will lead to different measurements for both normal and lateral forces as well as their fluctuations. 
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Appendix A: Direct calculation of lateral force fluctuations for two slabs with €2 — ei = e m 

Let us consider the lateral force fluctuations in the slab system in the absence of dielectric discontinuities, i.e. when 
e 2 = e i — e m- The standard three dimensional Coulomb interaction can be written as 

E=U*ufy r , (Al) 

2 J 47re e m [(x-y)2 + /2]5 

and the change in energy is thus 

SB = ( dxdy a -VPiW^(y) r , (A2) 
J 4^ e m [(x-y)2 + /2]2 

The average value of 6E is clearly zero but one can show that for delta-correlated charge distributions 

/XT? 2\ 9is92s f , , (a - (x-y)) 2 

From this one can extract the force correlator as 

where the integral over z is over the relative position x — y and the leading order term is proportional to A (there will 
be a correction term proportional to the perimeter dA of the region containing the charge on plate 1). The integral 
in Eq. (|A4[) is then easily evaluated to recover the result Eq. (|25p. 
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